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Abstract 

We investigate single top production in the presence of anomalous Wtb cou- 
plings. We explicitly show that, if these couplings arise from gauge invariant 
effective operators, the only relevant couplings for single top production and de- 
cay are the usual 7 M and a^ u q u terms, where q is the W boson momentum. This 
happens even in the single top production processes where the Wtb interaction 
involves off-shell top and/or bottom quarks. With this parameterisation for the 
Wtb vertex, we obtain expressions for the dependence on anomalous couplings of 
the single top cross sections, for (i) the t-channel process, performing a matching 
between tj and tbj production, where j is a light jet; (ii) s-channel tb production; 
(iii) associated tW~ production, including the correction from tW~b. We use 
these expressions to estimate, with a fast detector simulation, the simultaneous 
limits which the measurement of single top cross sections at LHC will set on V t b 
and possible anomalous couplings. Finally, a combination with top decay asym- 
metries and angular distributions is performed, showing how the limits can be 
improved when the latter are included in a global fit to Wtb couplings. 



1 Introduction 



In addition to top pair production, which is the largest source of top quarks in hadron 
collisions, single top processes will be of great importance for the study of the top quark 
properties at the Large Hadron Collider (LHC) [1]. With a cross section about three 
times smaller than for ti, single top quarks will also be coupiously produced through 
the electroweak Wtb interaction and, precisely because of this production mechanism, 
single top processes will provide essential information about the Wtb vertex. Their 
cross sections are proportional to the size of this interaction, and thus from their 
measurement the value of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element V t b 
will be obtained, as well as bounds on possible anomalous couplings [2,3]. Moreover, 
single top quarks will be produced with a high degree of polarisation [4] (in contrast 
to QCD ti production), which allows to build top spin asymmetries of order 0.1 — 0.4 
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in the leptonic decay t — > Wb — » ivb. Such asymmetries complement the observables 
measured in the W rest frame [5] in order to test anomalous Wtb couplings with a high 
precision and in a model-independent way. 

New physics beyond the Standard Model (SM) is expected to affect especially the 
top quark, and, in particular, it may modify its charged current interaction with its 
SU(2)l partner the bottom quark. For on-shell t, b and W, the most general Wtb 
vertex involving terms up to dimension five can be written as [5] 

£™ b = --^=bY(V L P L + V R P R )tW- 

-jf ^ (g L P L + g R P R ) t W~ + H.c. , (1) 

with q = p t — pb (being p t and pb the momenta of the top and b quark, respectively, 
following the fermion flow), which equals the W boson momentum. Additional a^ v k v 
and k^ terms, where k = p t + p b , can be absorbed into this Lagrangian using Gordon 
identities. If the W boson is on its mass shell or it couples to massless external fermions 
we have q^e^ = 0, where e M is the polarisation vector of the W boson, so that terms 
proportional to q^ can be dropped from the effective vertex. Within the SM, the 
only Wtb interaction term at the tree level is given by the left-handed 7 M term, with 
Vl = Vtb — 1. The rest of couplings are called "anomalous" and vanish at the tree level, 
although they can be generated by radiative corrections. They are not necessarily 
constants but rather "form factors", usually approximated by the constant term (as we 
will do in this work). If we assume that CP is conserved in the Wtb interaction then 
Vl,r and ql,r are real, and Vl can be taken to be positive without loss of generality. 

For off-shell top and/or bottom quarks the Lagrangian in Eq. (00) is not the most 
general one, and in principle it should be extended with k^ and a^ u k v terms. Nev- 
ertheless, if Wtb anomalous couplings arise from gauge invariant effective operators, 
single top production and decay can be described in full generality using the on-shell 
Lagrangian in Eq. (pQ) for the Wtb vertex, even in the processes where the top and 
bottom quarks involved in the Wtb interaction are far from their mass shell. This 
is a particular case of a more general set up, and it will be explicitly proved for the 
specific case of single top production. In particular, we will show that new physics 
contributions to the Wtb vertex can be "rewritten" using Gordon identities into the 
form of Eq. (pQ), even for t and b off-shell. The precise meaning of this "rewriting" 
will be clear in the next section: using an adequate parameterisation for the most 
general Wtb vertex, including k^ and a^ky terms, we will see that off-shell effects of 
Wtb anomalous couplings in single top production identically cancel when summed to 
contributions from anomalous gWtb quartic couplings, which are related to the former 
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by gauge invariance. We refer the reader to section [2], where these issues are explained 
in detail. 

After showing that for single top production the general Wtb vertex can be ade- 
quately described with the on-shell Lagrangian in Eq. (pQ), we will obtain expressions for 
the single top cross sections in terms of the couplings in Eq. (ED). There are three single 
top (and antitop) production processes in hadron collisions: (i) the t-channel process, 
also denoted as Wg fusion, which involves tj and tbj production; (ii) tb production 
with a s-channel W; (iii) associated tW~ production, with a correction from the tW~b 
process. In the SM their cross sections are proportional to V^(= Vj). For a general 
Wtb interaction the cross sections include additional terms involving anomalous Wtb 
couplings and can be conveniently written, factorising the SM cross section (calculated 
for Vl = 1), as 

<J = ^sm {Vl + ^ Vl + V l V r + g\ + g% + K ™* g L g R + ...), (2) 

where the k factors (9 numbers for the different combinations of couplings in £$| 6 ) 
determine the dependence on anomalous couplings. These factors are in general differ- 
ent for t and t production, and must be evaluated with a numerical integration of the 
corresponding cross section. They depend on parton distribution functions (PDFs), 
the factorisation scale Q 2 and parameters such as the top and bottom quark masses. 
These dependences translate into theoretical uncertainties which have to be considered 
when deriving limits on anomalous couplings from single top cross sections. 

We will then take the next step and estimate the potential limits on Wtb couplings 
from single top cross section measurements at LHC. These limits are not very stringent 
due to the experimental errors on the cross sections, which are expected to range 
between 13% (for the t-channel process) and 21% (for tb production). Finally, we will 
perform a combination with top decay observables, such as angular distributions and 
asymmetries, which complement the former and allow us to improve limits on Wtb 
couplings significantly, obtaining at the same time a measurement of V t b and limits on 
anomalous couplings. 

The rest of this paper is organised as follows. In section [2] we derive the Wtb 
couplings from gauge invariant effective operators and discuss the cancellation among 
off-shell contributions involving triple and quartic vertices. Readers mainly interested 
in the numerical results may skip this section. In sections EMS we study in turn the 
three single top production processes in the presence of anomalous couplings and de- 
rive the corresponding k factors for single top and antitop production. In section [6] we 
present a first estimate for the limits on anomalous couplings which can be obtained 
from the measurement of single top cross sections, as well as with their combination 
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with top decay angular distributions and asymmetries. Our conclusions are presented 
in section [7J The explicit proof of the cancellation among off-shell and quartic contri- 
butions is given in appendix [A] A brief review of the top decay observables used in 
the combination is given in appendix [Bj The effect of the top decay is examined in 
appendix O 



2 Wtb anomalous couplings from effective operators 

We follow the notation of Ref. [8] for gauge invariant effective operators with slight 
normalisation changes and sign differences. For reference, we summarise here the defi- 
nitions needed in this section. We denote by 

<?l = ( \ L ) , t R , b R (3) 



b 



L 



the weak interaction eigenstates for the third generation. The covariant derivative is 

D, = d, + ig^Gl + ig^Wl + ig'YB, , (4) 

where A a are the Gell-Mann matrices with a = 1 . . .8, t 1 the Pauli matrices for I = 
1, 2, 3, Y is the hypercharge and G^, and are the gauge fields for SU(3), SU(2)t 
and U(l)y respectively. The charged W boson fields are 

W ± = — ^= (W 1 =F iW 2 ) (5) 
V 2 

and the field strength tensor for SU(2)l is 

K» = W - d »K - 9^ijkW;W u k . (6) 
The SM Higgs doublet (ft has vacuum expectation value 

«>-*(") ' 

with v = 246 GeV, and we define <p = ir 2 (p*. 

The anomalous terms involving V R , g R and gi in the on-shell Lagrangian of Eq. ([T]) 
can arise from the effective operators 



! ' r2 b Rl n R w u 



t 



H = i{^iT 2 D^){i R ^b R ) - 

O uW = ((Il^tHr) 4> Wl u -> 2vb L %a^q v t R W~ , 

O dW = (qL^r'dn) W* u - [-2vb R ia^q v t L W~] f . (8) 
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Contributions to Vj, can originate from 

Of q = i^D^cf)) (ft/yMfe) 



<JV b Ll H L W u 



V2 



(9) 



Notice that these four operators do not yield quartic terms relevant for the amplitudes 
involved in single top production. Contributions to the Wtb vertex involving the top 
and b quark momenta can be obtained from the operators 

Du ={q L D tl t R ) D^, 



Du ={D^q L t R ) D»$, 
Dd = (q L D ll b R )D^, 
Dd = (D^q L b R )D^. 



(10) 



The opposite sign combinations 
Odu — Odu — * 
Ood — 0[) d — ► 



y b L kH R W~ - gg.v 5 L y g^t R W~ G" 
y b R kfH L W~ - gg s v b R ^t L W~ G" 



it 



(11) 



give k^ terms plus quartic gWtb interactions, while the same sign combinations 



^-b Lq n R w-, 



- 9 ^rb R qn L W u 



(12) 



give g M terms in the Wtb vertex but not gWtb ones, which cancel in the sums. Analo- 
gously, from the operators 

0' Du = i{q L e^D v t R ) Dp (f> , 
0' Gu = i(D v q L a^t R ) D,4>, 

0> Dd = i{D 1 ,q L <T> a 'b R )D ll <l>, (13) 

which are equivalent to the ones in Eq. ([TO for t, b on-shell [8], we can get contributions 
to terms. The combinations 



o' Du - 0' Du 



O' - O'- 

w Dd w Dd 



^ b L ia^Kt R W~ - gg s v b L ^ia^t R W~ G% , 



b R ia^k u t L W~ + gg s v bJ^ia^t L W~ G a v 



(14) 
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yield a^ u k v terms in the Wtb vertex plus quartic gWtb interactions, and the combina- 
tions 



0' Du + 0' Du - 9 -^-b L ia^q v t R W- , 



0' Dd + 0' Dd 



gv 



»—b R ia^q v t L W; (15) 



t 



give a^qu terms (which can also be obtained from the operators O u w and Odw m 
Eq. (jHJ), respectively) but not quartic gWtb ones. 

The most general Wtb vertex for t, b off-shell and W on-shell or coupling to external 
masless fermions includes the Lagrangian in Eq. (pp) plus o^ v k v and W terms. But, 
instead of simply adding the latter terms to the on-shell Lagrangian, it is much more 
convenient to parameterise the general Wtb vertex in terms of Cy^ tb plus operators 
which vanish when t and b are on their mass shell, 



Cwtb = Cwtb + Oi + 2 , (16) 



with 



O x = --J—b[ta^K(f 1L P L + f 1R P R )-(m b f 1L ~m t f 1R )YPL 

-{-mtfiL + m b f 1R )YPR ~ ^UxlPl + HrPr)] t W~ + H.c. , 

2 = --J—b[k^f 2L P L + f 2R P R )-ta^q,(f 2L P L + f 2R P R ) 
V 2,M W 

-(m b f 2L + m t f 2R )^P L - (m t f 2L + m b f 2R )^P R ] t W~ + H.c. , (17) 

being fiL,m, fiL,2R arbitrary constants. We point out that 0\, 2 are just the Gordon 
identities for the a^ v k v and k^ terms properly normalised. This parameterisation is 
completely general for a CP-conserving Wtb vertex. From Eqs. (fTTi l. (fT4l ) we observe 
that corresponding to the off-shell operators 0i >2 , the anomalous gWtb interactions 

(4) = (/^ + t htg; + h.c. , 



M w 2 



Of = ^f^b'-g^ {f 2L P L + f 2R P R )tW^G a u + Wc. (18) 

have to be introduced if the former arise from gauge invariant effective operators. 

A very important consequence of gauge invariance is the cancellation among the 
contributions of 0\ %2 and Ofl to the amplitudes. For the case of gb — > tW~ this is 
simbolically depicted in Fig. [D We emphasise that the cancellation takes place when 
all the terms in 0\ (or 2 ) written in Eqs. (fTTj) are included, not only the a^ u k v 
(or k^) ones. This clarification is important because only the a^ u k v and k^ terms in 



V2gg s T X 
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b W b W b w 




= 



9 t 9 t 9 t 



Figure 1: Cancellation between contributions to the gb — > tW amplitude. The gray 
circles in the diagrams stand for: (i) a Wtb interaction through 0\ in the first and 
second diagrams, and a gWtb interaction through in the third one; or (ii) the 
same involving 2 and 0% . 

have associated quartic vertices in Eq. (JTBl) . but choosing the parameterisation in 
Eq. ( |T7l ) makes the cancellation apparent, decreasing the number of couplings relevant 
for single top production from 8 to 4. The on-shell Lagrangian in Eq. can then be 
used in full generality for the study of tW~ production [9], as well as in the tbj and 
tW~b processes, for which the Wtb vertex involves off-shell t and/or b but the effects 
of 01,2 and Of\ cancel. (For tj and tb the top and bottom quarks are already on- 
shell.) The explicit proof of the cancellations is given in appendix [A] for the gb — > tW~ 
and gg — > tW~b processes. It is interesting to point out that in the latter case it 
implies that off-shell operators do not influence the top decay, and for ti production 
with decay of both quarks, i.e. gg — ► W + bW~b, the cancellation is expected to take 
place as well, although it is not proved explicitly. For tbj the diagrams can be related 
to tW~ production by crossing and the addition of a fermion line to the external W 
boson, and the proof is the same. 

A consequence of the gauge cancellation is that, as we have mentioned in the in- 
troduction, we can actually use Gordon identities to rewrite new physics contributions 
into the form of Eq. ((TJ), even for t and b off-shell. Let us assume that some new physics 
at a high scale gives a contribution to the Wtb vertex of the form 

AC = - ^ b W(c L P L + c R P R ) t W~ + H.c. , (19) 

with cl, Cr constants, plus additional triple and quartic contributions implied by gauge 
invariance. We can also write AC as 

AC = Ad + AC 2 , (20) 
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with 

A£i = -^b[W(c L P L + c R P R )-ia^q u (c L P L + c R P R ) 

-(m b c L + m t c R )^P L - (m t c L + m b c R )^P R ] t W~ + H.c. , 

A£ 2 = ~ b [i(J^q v {c L P L + c R P R ) + {m b c L + m t c R )^P L 

+ (m t c L + m b c R )YP R ] t W~ + H.c. (21) 

Because of the cancellation mentioned, the contribution of ACi (which has the same 
form as C 2 , with gf\L,ih = (2M w v / A 2 ) cl jR ) to single top production and decay van- 
ishes, thus we can effectively make the replacement 

AC -> A£ 2 , (22) 

which amounts to using the Gordon identities on AC independently of whether the 
top and bottom quarks are on-shell or not. Obviously, the same statement applies to 
cr^fcy-type contributions as well. 



3 The ^-channel process 

This process involves the scattering of a light quark and a gluon from the proton sea, 
as shown in Fig. [2] (a,b). (We will often refer to the tbj and tj processes generically, 



U d u d u d 




a b t g t b b t 



(a) (b) (c) 

Figure 2: Sample Feynman diagrams for single top production in the t-channel process. 
Additional diagrams are obtained by crossing the light quark fermion line, and/or 
replacing (it, d) by (c, s). The diagrams for antitop production are the charge conjugate 
ones. 

understanding that the same applies to the charge conjugate processes tbj and tj. 
When there is some difference we will note it explicitly.) Diagram (a) where the gluon 
splits into a bb pair is the dominant one in the SM, but diagram (b) has to be included 
as well in order to form a gauge invariant set. For low transverse momentum p\ of the 
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b quark, there are large logarithmic corrections to the cross section from the internal b 
propagator, which can be resummed by introducing a b quark PDF in the proton [10] 
and describing the process as 2 — > 2 body, as shown in diagram (c). The extra b 
quark is generated in this case from initial state radiation (ISR) in a parton shower 
Monte Carlo. When merging the 2 — > 2 and 2 — > 3 processes, one has to be careful to 
avoid double counting in the low p b T region. One possibility is to perform a matching as 
suggested in Ref. [11], based on the p b T distribution shown in Fig. [3] For the tbj process 
it is required that the b quark has transverse momentum above certain value, p h T > p™*, 
while for the tj process a veto p^ SR < p™* is imposed in the parton shower Monte Carlo. 
The cross section for the low-p^ region is determined from the next-to-leading order 
(NLO) value and the cross section above the cut, for which the perturbative calculation 
is reliable. This amounts to normalise the cross section in the low-p^ region with a K 
factor, 

Ka(tj, p^ SR < p c ^) = a NLO - a(tbj, p\ > p c ^) . (23) 

As it has been shown in Ref. [12], the choice p™ 1 = 10 GeV leads to a smooth transition 
between the high-p^ region from diagrams (a,b) and the low-p^ region, described by 
diagram (c) plus an additional &isr quark from ISR. Furthermore, it gives very good 
agreement for the kinematical distributions of various variables with the full QCD NLO 
calculation, as implemented in the generators ZTOP [13] and MCFM [14]. This is important 
in order to obtain an accurate prediction for the n coefficients corresponding to a^ u 
couplings, in which the vertex involves a momentum factor. This matching prescription 
also has the advantage that it avoids introducing negative weight events, as it happens 
with procedures involving a subtraction term [15, 16] and can easily be implemented 
in a Monte Carlo generator. Electroweak corrections are small, only a —1.5% [17], and 
do not modify the kinematics. 

For each tj, tbj sub-process bu — > td, gu — > tbd, etc. (some of the amplitudes are 
related by crossing symmetry) the squared matrix elements are calculated analytically 
with FORM [18], taking the Wtb interaction in Eq. (PQ) and keeping nib non- vanishing. 
Then, the contributions to the cross section of the different products of anomalous 
couplings in Eq. j2j) are obtained by identifying in the squared matrix element the term 
multiplying each combination of couplings, and then integrating these terms separately 
over phase space and PDFs. Obviously, the cross section is the sum of all terms, and 
the k coefficients in Eq. ([2]) are obtained by dividing by the Vl term. A cross-check 
of all matrix elements is performed calculating them numerically with HELAS [19], 
extended to include the non-renormalisable a^ v and couplings in Eqs. (pQ) and (TT71) . 
Both the analytical and numerical matrix elements agree within several digits. For the 
calculation of the cross sections we use CTEQ6M [20] PDFs with the scales advocated in 
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10 ' 10 ' ' 20 ' ^ 30 ' 40 ' ' 50 

b 

Pt 

Figure 3: Normalised transverse momentum distribution for the b quark in the t-channel 
process, when it is generated from the hard process or from ISR in the parton-shower 
Monte Carlo. 

Ref. [21]: Q 2 = q 2 for the light quark, Q 2 = q 2 +m 2 for the initial b and Q 2 = {p b T ) 2 +ml 
for the gluon (q is the momentum of the t-channel W boson). We take m t = 175 GeV, 
M w = 80.39 GeV, m b = 4.8 GeV, and use the NLO cross sections cr NLO (t) = 155. 91U 
pb, ctnlo(^) = 90.7^3 7 pb [13] (the uncertainties have been rescaled to Am t = 1.8 
GeV). 

The cut p h T > 10 GeV in the tbj process can readily be performed at the generator 
level, obtaining cross sections cr{tbj) = 80.9 pb, critbj) = 48.3 pb. The veto p^ SR < 10 
GeV for the tj process is more involved. It is implemented by linking Pythia 6.4 [22] 
to the Monte Carlo integration program. During the numerical integration of the cross 
section, each weighted event is feeded into Pythia, which adds the ISR. Subsequently, 
the event is accepted (in which case it is considered for the integration) or rejected. 
The effect of the p^ SR < 10 veto on the k coefficients amounts to a few percent, thus 
the systematic uncertainty originating from ISR modeling in Pythia is expected to be 
small. The results for the most relevant coefficients in tj, ij, tbj and tbj are shown 
in Tables [Tj and [2J omitting for brevity the coefficients which are smaller than 0.1 for 
both t and i production. In each block, the first column corresponds to the interval of 
variation when CTEQ6M or MRST 2006 [23] PDFs are used. The second column is 
the uncertainty from the choice of factorisacion scale, and the third and fourth columns 
represent the uncertainties corresponding to Am t = 1.8 GeV and Amj = 220 MeV, 
respectively. Monte Carlo statistical uncertainties are of order 10~ 3 in most cases, 
and the same random seeds are used in all samples in order to reduce the statistical 
fluctuations among them. The most salient features of the results in Tables [T] and [2] 
are: 
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tj 










tj 








K 


AQ 


Amt 




K 




AQ 


Amt 


A?71b 




0.916-0.923 


+0. 

-0. 


+0. 
-0. 


+0. 
-0. 


1.082 - 


1.084 


+0. 
-0. 


+0. 
-0. 


+0. 
-0. 


A 


1.75 - 1.79 


+0.044 
-0.038 


+0.007 
-0.035 


+0 

-0.027 


2.16- 


2.17 


+0.035 
-0.022 


+0.014 
-0.032 


+0. 
-0. 


9l 


2.18 


+0.042 
-0.033 


+0.014 
-0.034 


+0. 

-0.022 


1.75- 


1.77 


+0.042 
-0.033 


+0.007 
-0.033 


+0. 

-0.025 


Vl9r. 


-(0.348-0.365) 


+0.007 
-0.011 


+0. 
-0. 


+0. 
-0. 


-(0.038- 


- 0.040) 


+0.010 
-0.009 


+0. 
-0. 


+0. 
-0. 


V R g L 


-(0.006-0.008) 


+0.006 
-0.005 


+0. 
-0. 


+0. 
-0. 


-(0.399- 


- 0.408) 


+0. 

-0.008 


+0. 
-0. 


+0. 
-0. 



Table 1: Representative k factors for the tj and tj processes and their uncertainties, 
explained in the text. Errors smaller than 0.005 are omitted. 







tbj 








tbj 








K 


AQ 


Am t 


Amt 


K 


AQ 


Am t 


A 771 (, 


V 2 


0.927- 0.932 


+0.005 
-0. 


+0. 
-0. 


+0. 
-0. 


1.068- 1.069 


+0. 

-0.005 


+0. 
-0. 


+0. 
-0. 


V L V R 


-0.117 


+ 0. 
-0. 


+0. 
-0. 


+0.005 
-0.005 


-0.126 


+0. 
-0. 


+0. 
-0. 


+0.006 
-0.006 


al 


1.96-2.01 


+0.070 
-0.056 


+0.005 
-0.005 


+0. 
-0. 


2.98-3.00 


+0.040 
-0.040 


+0.014 
-0.014 


+0. 
-0. 


9k 


2.97-2.98 


+0.056 
-0.043 


+0.013 
-0.013 


+0. 
-0. 


2.08-2.11 


+0.056 
-0.045 


+0.006 
-0.007 


+0. 
-0. 


Vl9r 


-(0.539- 0.550) 


+0.012 
-0.010 


+0. 
-0. 


+0. 
-0. 


-(0.169 - 0.172) 


+0.010 
-0.010 


+0.014 
-0.013 


+0. 

-0. 


V R g L 


-(0.121 - 0.134) 


+0.009 
-0.011 


+ 0. 

-0. 


+0. 
-0. 


-(0.567- 0.571) 


+0.014 
-0.013 


+0. 

-0. 


+0. 
-0. 



Table 2: Representative k factors for the tbj and tbj processes and their uncertainties, 
explained in the text. Errors smaller than 0.005 are omitted. 

(i) The coefficients of the terms are not equal to unity, in contrast with the other 
single top processes studied in the next two sections. This is easy to understand 
with the examination of the squared matrix element for the bu —>■ td and ud — ► tb 
processes, related by crossing symmetry. In both cases the coefficient of V£ 
is proportional to (p u ■ Pb)(Pd • Pt) (see Fig. EJ), while the coefficient of is 
proportional to (p u ■ Pt){pd • Pb)- In the s-channel process ud — > tb, integration in 
phase space renders the coefficients of V£ and equal, but this is not the case 
for the t-channel process bu — > td. In the latter, the two coefficients are different 
even in the case of all fermions massless. 

(ii) The coefficient of the term is different for single top and single antitop pro- 
duction, but the differences cancel to a large extent in the total cross section. 
This property makes the ratio R(t/t) = a(t)/a(t) more sensitive to a Vr compo- 
nent than the total cross section itself. A purely left-handed interaction yields 
a total (top plus antitop) cross-section of 246 pb, while a purely right-handed 
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interaction gives a total of 241 pb. Even in this extreme case, the difference 
is only of 1.9%, too small to be observed given the experimental uncertainties 
involved (see section [6]). However, for a left-handed interaction R(i/t) = 0.58, 
while for a right-handed interaction R(i/t) = 0.68, and this 17% difference should 
be visible. Nevertheless, its sensitivity to Vr is superseeded by top decay angular 
distributions [6,7] (see appendix IBl) I 1 ! 

(iii) The coefficients of the g\ and g R terms are larger than unity, as expected from the 
q v enhancement factor in the interaction. Moreover, they are larger for tbj than 
for tj, due to the larger energy involved when p 1 ^ is larger (hence the importance 
of describing correctly the tail of the "p h T distribution). 

(iv) Interferences among anomalous couplings are important, and in some cases the 
corresponding terms have coefficients of order unity. This implies that taking only 
one nonzero anomalous coupling at a time is by far a too simplistic assumption, 
and the possible cancellations among anomalous couplings have to be explored 
and constrained using as much information as possible, from single top production 
cross sections as well as from top decay asymmetries. 

We conclude this section by comparing these results with the ones obtained using 
an alternative procedure to remove double counting, which is to subtract from the 
cross section a term involving the first order logarithmic corrections included in the 
2^3 process [15,16]. This term <r su b is calculated by replacing in the calculation of 
the bq — > tj cross sections the b quark PDF by 

f g (x/z,Q), (24) 

where f g is the gluon PDF. For convenience, the double counting term can be subtracted 
from the tj contribution. The cross sections obtained can be read in Table [3] In order 
to compare with the results obtained using the matching prescription we separate the tj 
and tbj cross sections in two regions, depending on whether -p h T is larger or smaller than 
10 GeV. The total cross section is larger than in previous case, where we normalised 
cross sections to the NLO value. Notice that with the subtraction method a 20% of 
the cross section for p b T > 10 GeV is given by the 2 — > 2 process tj plus ISR, which has 
a softer spectrum. The k factors for each region can be calculated by summing the tj 
(minus the subtraction term) and tbj contributions for each coupling and dividing the 

1 A Vr term is indirectly constrained by the measured rate of b — > s-f [24-27] but, if our aim is 
to obtain direct and model-independent measurements of the Wtb interaction we cannot rely on such 
indirect measurements which, on the other hand, are much more restrictive. 



~ a s {Q) Q 2 f'dz 



12 



result by the V£ term. The most important k factors are listed in Table [H calculated 
with this method and with the matching at 10 GeV. We also include the values obtained 
using only the tbj process, for both the low and high p b T regions. 

p\ < 10 GeV v\ > 10 GeV Total 

(r(tj)-a enb (tj) 28.3 16.4 44.6 

a{tbj) 50.4 80.9 131.3 

v(tj) - a mh (tj) 18.1 10.8 28.9 

a(tbj) 31.6 48.3 79.9 

Table 3: Cross sections (in pb) for the t-channel process when the subtraction procedure 
is applied. 
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S tbj 
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S tbj 




0.915 


0.921 0.916 


0.927 


0.923 0.927 


9l 


1.75 


1.61 1.60 


1.96 


1.98 1.96 


9n 


2.18 


1.99 2.01 


2.97 


2.90 2.97 



Single antitop 





Pt 


< 10 GeV 


PT 


> 10 GeV 




M 


S tbj 


M 


S tbj 




1.084 


1.075 1.081 


1.068 


1.073 1.068 


9l 


2.16 


2.03 2.11 


2.98 


2.90 2.98 


9l 


1.75 


1.66 1.71 


2.08 


2.08 2.08 



Table 4: Comparison of some n factors calculated with the matching (M) and subtrac- 
tion (S) prescriptions, and using only the tbj process 

For k Vr the agreement is very good, better than 1%. For k 9l and k 9r the values 
obtained with the subtraction mathod are generally smaller. The agreement for p T > 10 
GeV, where the dependence on anomalous couplings is stronger, is better than 3%, of 
the same order of the theoretical uncertainties quoted in Table El The differences for 
p T < 10 GeV are more significant but always below 10%. In summary, both methods 
for removing double counting give results not very different, and similar to the ones 
obtained just using the tbj process. The results obtained with the matching are believed 
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to be more accurate, for the reasons explained above. 



4 tb production 

This process is mediated by s-channel Feynman diagrams as the one depicted in Fig. [H 
In Ref. [13] it has been shown that NLO corrections do not significantly modify the 
kinematics and can be accounted for by a K factor. Thus, the effect of anomalous 
Wtb couplings is expected to be well approximated by including this K factor in the 
tree-level cross section involving anomalous couplings. The NLO cross sections are 
Cnlo(^) — 6.56±0.38 pb, cx NL o(t) = 4.09±0.24 pb [13]. Partial higher-order results are 



u b 




d t 



Figure 4: Feynman diagram for single top production in the tb process. An additional 
diagram is obtained replacing (u, d) by (c, s). For antitop production the diagrams are 
the charge conjugate ones. 

also known [28], and increase the tb cross section by about ~ 10%. Their impact in the 
tb final state experimentally analysed, consisting of tb production plus contamination 
from tj and tbj (see section [6]) is of only 3%, much smaller than the experimental 
uncertainty which is around 20%. Anyway, they can be included with a K factor. 

In our calculations we use CTEQ6M PDFs with Q = y/s, obtaining a(tb) = 4.62 
pb, a{tb) = 2.88 pb. The calculation of the k factors and their uncertainties proceeds 
in the same way as for tj and tbj in the previous section. Results are shown in Table El 
including uncertainties from PDFs, factorisation scale, the top and the b quark masses. 
We point out that the factors multiplying V£ and in the cross section are equal, 
then k Vr equals unity. We also have k 9l = k 9r , k Vl9l = k Vr9r , k Vl9r = k Vr9l . From 
Table E] we observe that: 

(i) The k factors of g\ and g\ are a factor of four larger than for the t-channel process, 
because in tb production the s-channel W boson carries a larger momentum, and 
so the q v factor in the a^ u vertex gives a larger enhancement. 
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-0. 


12.22-12.28 


+0.21 
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-0.008 
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-0.012 
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1.25 


+0.008 
-0.009 


+0.013 
-0.013 


+0.056 
-0.056 


Vl9l, V R g R 


-0.415 


+0. 

-0. 


+ 0. 

-0. 


+0.018 
-0.018 


-0.426 


+0. 
-0. 


+ 0. 

-0. 


+0.019 
-0.019 


Vl9r, V R g L 


-5.51 


+0.009 
-0.010 


+0.057 
-0.057 


+0. 
-0. 


-5.48 


+0.008 
-0.010 


+0.057 
-0.056 


+0. 
-0. 



Table 5: k factors for the tb and tb processes and their uncertainties, explained in the 
text. Errors smaller than 0.005 are omitted. 

(ii) For tb and tb production the factors are very similar, although not equal (the 
difference is not due to Monte Carlo statistics, which is very high). Then, the 
measurement of the ratio a(tb)/a{tb) is not as useful as in the t-channel process. 

(iii) Interferences among couplings are again important, in particular between Vl and 
gu, and between Vr and g^. 

Finally, it must be remarked that tb production, with a cross section much smaller than 
the t-channel process, has a large contamination from the latter in a real experiment. 
This is taken into account in the limits presented in section [6l 

5 tW production 

At lowest order, the gb — > tW~ process is mediated by the two diagrams in Fig. El where 
the initial b quark comes from splitting g — > bb of a gluon in the proton sea. The charge 
conjugate process exhibits the same features and has the same cross section, so we will 
often refer only to tW~ production for brevity. Even at this order the t and b quarks 

b W b W 




9 t 9 



Figure 5: Feynman diagrams for single top production in the gb — > tW process. 
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entering the Wtb vertex are not on their mass shell in the first and second diagrams, 
respectively. Like the t-channel process, tW~ production has a NLO correction from 
gg — > tW~b. However, an important difference is that the full gauge invariant set of 



Figure 6: Feynman diagrams for gg — > tW~b. The two initial gluons are labelled as g\ 
and #2 for clarity. The last three diagrams correspond to resonant ti production with 
decay i — > W~b. 

8 Feynman diagrams for gg — > tW~b (see Fig. EJ) also includes on-shell ti production 
with t — > W~b, whose cross section is about one order of magnitude larger than for 
tW~ itself. Apart from the displeasing fact of considering ti as a huge "correction" to 
tW~ production, this is not convenient from a practical point of view. The motivation 
for studying single top production in this paper, and perhaps the main reason to study 
these processes at all, is to measure V tb and set limits on anomalous Wtb couplings. 
But the ti cross section is actually independent of the Wtb interaction. Top pairs are 
produced in a QCD process and, as long as t — > Wb dominates the top decay (which 
is expected, due to the smallness of V t d and V ts ), the branching ratio is independent as 
well. It must be pointed out, however, that if one restricts the invariant mass of the 
W~ boson and b quark m-wb to select only off-peak contributions, for example requiring 
\m Wb — m t \ > 20 GeV, then the gg — > ti — > tW~b cross section in this region does 
depend on the Wtb coupling, e.g. it scales with Vl like in single top production. Still, 
the total cross section in all phase space remains independent, because the cross section 
around the peak (whose height is determined by the top width r t , which depends on 
Wtb couplings) compensates the variations in the region where the i quark is off-shell. 

Since the total gg — > ti —> tW~b cross section is practically independent of the 
Wtb interactional if one wants to study the influence of the latter on the total tW~b 
cross section it is natural to consider gg — » ti — >• tW~b as "background" and the excess 

2 Of course, this does not preclude the fact that the i — » W~b angular distribution and the W 
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of events over it, which depends on Vl and anomalous couplings, as "signal". This 
approach is convenient from the point of view of the experimental analysis as well. 
The single top signal is then 

Signal = <r%(gg -> tW~b) - a 3 (gg -> tt -> tW'b) , (25) 

where the first term a 8 includes the full set of 8 diagrams in Fig. [6] and the second term 
<7 3 only involves the three resonant ti production diagrams with i decay (the last three 
diagrams). Clearly, 

Signal = <rs(g9 -> tW~l) + CT int , (26) 

so that the tW~b excess over the tt "background" is due to the 5 non-resonant diagrams 
in Fig. E] (05) plus their interference with tt production (this equation defines a-^). This 
interference is not negligible, and amounts to a —20%. Note that this calculation of 
the single top contribution involving subtraction of tt involves small violations of gauge 
invariance of order T t /m t . 

There is a different (gauge invariant) method in the literature [42,43] proposed to 
remove the tt contribution from tW~b, which consists in performing a cut requiring 
\mwb — m t\ larger than some quantity, of the order of 15 — 20 GeV. But still with this 
kinematical cut the tt contribution is comparable to the non-resonant one. In Table [6] 
we collect, for various m Wb cuts, the cross sections corresponding to several subsets 
of diagrams, using CTEQ6L1 PDFs with Q = m t + M w [44] § Notice that even for 
\mwb — fnt\ > 50 GeV the ti contribution 03 is comparable (a 40%) to the non-resonant 
one cr 5 , and their interference is large (—30%). This is not surprising, because the 
three ti diagrams do not have any particular suppression, and outside the mwb peak 
they are expected to have similar size as the non-resonant ones. Thus, ti contributions 
cannot be fully removed by this method (see also Ref. [45]). One could still argue 
that ti production outside the my/b peak behaves as single top production, because the 
cross section in this phase space region depends on Vl and anomalous couplings, and 
might be tempted to include it as single top production. But separating ti — > tW~b in 
two parts, near the peak and off-peak, makes the former (which is judiciously taken as 
"background") also dependent on Wtl i couplings and on, y complicates the analysis^ 



polarisation, as well as the kinematical distributions determined by them, depend on the structure of 
the Wtb vertex, as it is well known [3,5,29-41]. 

3 To our knowledge, there is no study in the literature comparing the leading order (LO) and NLO 
kinematics for tW~ production, unlike in s- and i-channel single top production [12,13]. Then, we 
conservatively use LO PDFs for the tW~ and tW~b processes. 

4 Note however that a cut on mwb will likely be useful to reduce the ti background and thus the 
experimental uncertainty associated to the subtraction. Then, in the real analysis a it subtraction 
with a mwb cut may yield the best results. This issue requires a detailed evaluation of systematic 
uncertainties and is beyond the scope of the present work. 
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no cut 439 426 

15 GeV 24.7 13.6 

25 GeV 18.0 8.20 

50 GeV 10.70 4.19 



08 - 03 


05 


0int 


12.9 


16.5 


-3.6 


11.1 


14.5 


-3.5 


9.81 


13.2 


-3.4 


6.51 


9.71 


-3.20 



Table 6: Cross sections (in pb) for gg — > tW~b involving different subsets of Feynman 
diagrams: the full set (<7 8 ), resonant ti production (cr 3 ), non-resonant diagrams (ct 5 ) 
and interference between ti and non-resonant diagrams (a- mt ). A kinematical cut on 
\mwb — mt | is applied in all but the first row (see the text). For the charge conjugate 
process gg — > iW + b the cross sections are equal. 

There is a very interesting feature which can also be observed in the last column of 
Table [6l The net interference between ti — > tW~b and non-resonant diagrams nearly 
vanishes around the mwb peak. Interference terms have the form 

ReMtM, = ReM*M 2 _ J, , (27) 

m wb m t + t 

where .Mi, .Mj are the amplitudes of a non-resonant and resonant diagram, respec- 
tively, and M. j the latter without the i propagator, which has a momentum p| = m^ 6 . 
Since the two amplitudes do not have relative complex phases (except for the i propa- 
gator) the product M.^ = A4*Aij is real, and we have 

i 2 2 

777 777 

Re M*Mj = MijRe — = M y -^ Wfc f . (28) 

mfo-mt + im t T t {m l Wh - mf) 2 + (m t T t ) 2 

Near the peak at mf the factor Jidy is expected to be approximately constant, so that 
integration in m\y h within a symmetric interval [m 2 — nT t , m 2 + nT t ] with n not too 
large gives a vanishing interference. The cancellation can be nicely observed in Fig. [71 
for mwb < frit the interference is positive, while for my/b > frit it is negative, leading 
to a non-trivial peak-dip structure. This also implies that a- mi is rather independent 
of the top width, even if it is changed by a factor of four§| This feature is crucial for 
our analysis, because it allows us to safely neglect the dependence of T t on anomalous 
couplings, and justifies expanding as — a 3 in terms of K factors and products of Wtb 
couplings, as done in Eq. (j2J). (Non-resonant contributions in 05 are already practically 
independent of r t , which is negligible in the denominators.) 

Subtracting the ti — > tW~b contribution and considering a 8 — a 3 as our signal also 
has some drawbacks. For example, one has to introduce negative weight events, as 



5 It must be noted that the effect of anomalous couplings in the top width is much more modest 
than in single top cross sections due to the different momenta scales involved. 
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Figure 7: Kinematical distribution of mwb for the single top signal in Eq. (1261) and its 
two contributions. The cross sections are in pb. 

it is apparent from Fig. [3 More importantly, the p b T distribution does not exhibit a 
good behaviour, and even becomes negative for > 80 GeV. In Fig. [8] we present 
this distribution for tW~b in three cases: (i) when the tt — > tW~b contribution is 
subtracted; (ii) when a cut of 15, 25 or 50 GeV is applied on \m Wb — m t \; (iii) when the 
b quark is generated from ISR added to the tW~ process by Pythia. The distribution 
after tt subtraction is very soft, in contrast with the case when a cut is applied on mwb, 
indicating that tt — > tW~b give the b quarks with higher even away from the peak. 

It 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 : 



— Subtraction 
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Figure 8: Normalised transverse momentum distribution for the b quark in the tW~b 
process in several cases: for the hard tW~b process subtracting the tt contribution 
(solid line), with a cut on \mwb — Trit\ of 15, 25 or 50 GeV (dotted lines, from top 
to bottom) and from ISR added to the tW~ process in a parton shower Monte Carlo 
(dashed line). 

This bad behaviour of the jP T distribution makes it difficult to perform a matching 



19 



between tW~ and tW~b like the one used in section [3] for the t-channel process, unless 
some unrealistic p™* is chosen. Instead, we remove double counting between tW~ 
and tW~b by subtracting the first order logarithmic corrections included in the 2^3 
process [43,46], calculated by replacing in the calculation of the gb — > tW~ cross section 
the b quark PDF by the function in Eq. ( |24l) . The subtraction term is rather large [43]: 
using CTEQ6L1 PDFs and Q = m t + M w we find a(gb -> tW~) = 28.5 pb, and for the 
subtraction term a su b = 21.0 pb. Notice that a(gb — ► tW~) — cr sah + cr 5 (gg — > tW~b) = 
24.0 pb, a 15% smaller than the first order result for a(gb — > tW~). The interference 
with tt — > tWb further decreases the total cross section to 20.4 pb, a 28% smaller than 
a(gb^tW-). 

For tW~ we calculate the k factors as in the previous single top processes, using 
a squared matrix element calculated analytically with FORM (and numerically checked 
with HELAS) in which the different contributions are identified. For tW~b, the larger 
number of diagrams (8, giving 36 terms in the squared amplitude) and the number 
of products of couplings in each make this approach unsuitable. We calculate the k 
factors numerically by evaluating the squared matrix element iteratively, with different 
Wtb couplings. Taking one coupling equal to one and the rest zero, we can obtain 
the four coefficients of quadratic terms. Then, taking two couplings equal to one 
and the rest vanishing allows to extract the six interference terms, subtracting the 
quadratic contributions previously calculated. The results are presented in Tables d 
and [8] omitting the k coefficients smaller than 0.1. The uncertainties have the same 
source as in the previous sections: the intervals in the first column correspond to the 
PDF uncertainty, and the last three columns to the uncertainties from factorisation 
scale, the top mass and the bottom mass, respectively. In these processes both top and 
antitop production have the same values of k, providing an additional cross-check of 
our calculations. 

k AQ Am t Amj 

Vi 1 

„2 _2 o / |( i _q C7 +0.23 +0.015 +0.009 

yLiyR °-^ u -0.11 -0.015 -0.008 

V L g R} V R g L 1 



Table 7: Representative « factors for the tW and tW + processes and their uncertain- 
ties, explained in the text. Errors smaller than 0.005 are omitted. 
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V L 9R,V R g L 1.21-1.23 



+0.19 
-0.04 
+0.014 
-0.003 



+0.009 
-0.027 



+0.005 
-0.007 



+0.030 



-0. 
+0. 
-0. 



Table 8: Representative k factors for the tW b and tW + b processes and their uncer- 
tainties, explained in the text. Errors smaller than 0.005 are omitted. 

6 Limits on anomalous couplings 

In this section we present estimates of the limits that can be obtained from single top 
cross section measurements, and also with their combination with top decay observ- 
ables, which can be measured either in single top or top pair production. We will show 
that significant bounds on anomalous couplings can be obtained despite the possibility 
of cancellations and the contamination among the three single top processes. In these 
estimates we will assume sensitivities for single top cross section measurements based 
on previous literature, and for top decay asymmetries A± and helicity ratios pr^l (see 
appendix [B]) we will take the experimental uncertainties from ti production. On the 
other hand, the experimental precision for the measurement of the asymmetry ratio r&; 
has not been estimated as yet. Nevertheless, the results are weakly dependent on the 
sensitivity of this observable, provided that it is better than ~ 8%. A more complete 
study will be presented elsewhere, including all experimental systematic uncertainties 
and the SM background. 

6.1 Limits from single top production 

The calculations of sections [3H5] have provided us with theoretical expressions for the 
single top cross sections including anomalous Wtb couplings. In order to link them 
to real or simulated experimental data, the detection efficiencies for each single top 
process (i.e. the fraction of events which survive the selection criteria required for 
the experimental analysis) must be known. First, because comparing an experimental 
sample with a theoretical cross section obviously requires to know which fraction of 
the events produced are actually present in the sample analysed. But also because 
different single top processes will contribute to a given final state, and the relative 
weights determine the dependence on anomalous couplings of the measured cross sec- 
tion. Efficiencies depend on the final state considered, selection criteria, etc. and must 
be determined at least with a fast detector simulation. Once these efficiencies are found 
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for each final state of interest, the experimental measurement of the cross section and 
its uncertainty allow to set limits on anomalous couplings. 

This is best explained with an example. Experimentally, the t-channel process can 
be investigated in final states with semileptonic top decay, requiring the presence of 
a charged lepton, significant missing energy, a 6-tagged jet (from the top decay), a 
forward jet with large transverse momentum and no additional jets with large pr to 
suppress tt production (see for example Ref. [47]). On the other hand, the theoretical 
calculation of this process is conveniently separated into tj and tbj production. Their 
respective efficiencies e t j and e t ij can be computed with a Monte Carlo simulation, 
including detector effects. Then, the predicted number of single top events is 

N = N tJ + N tlj = e tj a(tj) SM (vZ + KlfVi + K% v *V L V R + ...) 

+e tb] a(tbj) SM (Vl + «g Vl + Klf* V L V R + ...) , 
+e tb (T(fb) SM (Vl + V* + kY r V l V r + ...), 
+e tw a(tW-) SM (Vl + Vl + V l V r + ...), 
+e tWb a(tW~b) SM (Vl + Vl + V L V R + ...), (29) 

where we have introduced additional subscripts on the k factors to distinguish the 
different processes. The last three terms represent the "contamination" from the other 
single top processes, which can be reduced (i.e. making e t b, £tw, £twb "C e t j,e t bj) with 
suitable selection criteria. For antitop production the expression is the same, but 
substituting the numerical values for those corresponding to the charge conjugate pro- 
cesses. It is apparent that the functional dependence on Vl and anomalous couplings 
of the number of events observed is determined not only by the k factors but also by 
the SM cross sections and the efficiencies. This approach is different from the usual 
single top analyses: here all single top processes are considered as signals, with differ- 
ent efficiencies and dependence on anomalous couplings. Then, three final states can 
be selected trying to isolate each of the processes (t-channel, tb and tW~ /tW~b) as 
much as possible. Nevertheless, contributions from all processes will always be present, 
especially from the t-channel process which has a much larger cross section. 

A second issue to be considered is the dependence of efficiencies on Wtb couplings. 
Anomalous couplings change the angular and energy distributions both in the produc- 
tion and the decay of the top quark, and then they affect the efficiencies. A complete 
scan over the 3-dimensional space of anomalous couplings to parameterise the efficien- 
cies is involved, but a simpler approach can also be followed. First, the SM efficiencies 
can be used to obtain bounds on anomalous couplings. With these bounds the efficien- 
cies can be reevaluated if a positive signal beyond the SM is found and, in any case, 
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their uncertainties due to anomalous couplings must be estimated and included in the 
final results if they are found relevant. A similar procedure is foreseen for the study of 
anomalous couplings in top decays [7]. 

For these evaluations we have developed the Monte Carlo generator Protos for the 
different single top production processes considered, including the decay of the top 
quark and W boson(s) with finite width and spin effects. Matrix elements for 

bq -> tq' -> W + bq' -> l+vbq' , 
gq — > tq'b — > W + bq'b — > £ + ubq'b, 
qq' W + bb -> tvbb , 

gb -> tW~ -> W+ftW" -> £ + z/6gg' , 

-> *W~6 -> 6W6 -> tvbqqb (30) 

and their charge conjugate are calculated with HELAS, with VFtfr anomalous couplings 
implemented in the production as well as in the decay of the top quark. The output of 
the generators provide events with the colour information necessary for hadronisation 
in order to be interfaced to Pythia. These generators are also used to check the effect 
of the top quark decay in the k factors, which is discussed in appendix O 

We calculate the efficiencies e using a fast simulation of the ATLAS detector [48]. 
We restrict ourselves to the electron channel, final states with a muon are alike. Event 
samples for each the process are generated for final states containing at least one e ± , 
corresponding to luminosities of 30 fb _1 for tj and tbj, 300 fb" 1 for tb and 60 fb _1 
for tW~ and tW~b. For the latter process, this requires simulating about two times 
more events, so that the sum of events with positive and negative weight corresponds 
to the cross section times the luminosity. Results are finally rescaled to 30 fb _1 . A K 
factor of 1.54 [44,49] is introduced in the tW~ , tW~b processes for consistency with 
the other channels, resulting in cross section a(tW~ + tW~b) = 31.4 pb. Events are 
passed through Pythia including pile-up, ATLFAST and ATLFASTB, where a b tagging 
efficiency of 60% is chosen. We concentrate ourselves on a final state for each single 
top production process, requiring in all of them the presence of an isolated electron 
or positron with px > 25 GeV and missing energy larger than 25 GeV. The rest of 
selection criteria, adopted to isolate each process as far as possible, are respectively: 

(i) Final state 1 (for the t-channel process): a forward jet with pseudorapidity in the 
range 2.5 < \rj\ < 5 withp^ > 50 GeV; at least one central b jet wither > 30 GeV; 
at most one additional central jet, which cannot have px > 30 GeV. The top quark 
mass is reconstructed from the charged lepton and b quark momenta (if there are 
more than one b jet we select the one with largest transverse momentum), and 
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the missing energy. The transverse neutrino momentum is assumed to equal the 
missing energy of the event, and the longitudinal component is found solving a 
quadratic equation requiring that (p e + p u ) 2 = M^. The solution selected is the 
one with smaller p uz . Once that the top quark momentum and its invariant mass 
are reconstructed, the latter is required to be between 150 GeV and 225 GeV. 

(ii) Final state 2 (for tb production): we require that neither forward jets nor light 
central jets are present with transverse momentum larger than 15 GeV. Final 
states must have two 6-tagged jets with pt > 30 GeV. 

(iii) Final state 3 (for tW~ and tW~b production): the signal is enhanced requiring 
one (and only one) b jet with pt > 30 GeV and at least two light jets with px > 50 
GeV. The invariant mass of the W decaying hadronically is reconstructed from 
the two light jets with the largest transverse momentum, and it must be between 
60 and 110 GeV. 

The number of events for the five production processes considered are listed in Table [9], 
before and after the selection criteria for the final states 1 (tj), 2 (tb) and 3 (tjj), 
respectively. Comparing with Ref. [47], from where we take the experimental (statis- 
tical plus systematic) uncertainties expected, we observe that these selection criteria 
reproduce the features most important for our analysis: (i) the t-channel process, with 
a cross section much larger, can be cleanly separated from the other two processes; (ii) 
the tb process has a large contamination from i-channel single top production, which 
is about two times larger even after selection criteria are applied; (iii) tW~ + tW~b 
production can also be separated from the other processes, and is about 4 times larger 
after selection. Exact agreement with Ref. [47] is not expected, because the signal 





No sel 


sel 1 


sel 2 


sel 3 


tj 


391200 


7914 


398 


570 


tbj 


430800 


6231 


1411 


1265 


tb 


35500 


271.2 


957.7 


120.0 


tw- 


145200 


250.5 


106.0 


3455 


tW'b 


249800 


392.0 


111.5 


4869.5 



Table 9: Number of events for the different single top production processes in final 
states containing at least one e ± (first column) and after the selection criteria corre- 
sponding to the three final states studied. The luminosity assumed is 30 fb -1 . 

modeling is different (only tbj is considered in the t-channel process and the tW~b 



24 



correction is not added to tW~ production) and pile-up was not included in those sim- 
ulations. For these reasons we have not applied the same selection criteria but have 
slightly adapted them. Optimisation of the cuts is of course possible, as well as the use 
of a likelihood analysis to further suppress the SM background, mainly constituted by 
top pair production. This is beyond the scope of the present discussion. 

The expressions of the cross sections for arbitrary Wtb couplings have been imple- 
mented in TopFit [5], including their theoretical uncertainties. The efficiencies for the 
different final states are calculated from Table [9J and it is assumed that the number 
of measured events equals the SM expectation. The experimental errors assumed for 
cross section measurements in the three final states studied are [47]: 

Final state 1 (tj): 1.0% (stat) © 11% (exp) © 6% (bkg) © 5% (lum) , 
Final state 2 (tb): 12% (stat) © 12% (exp) © 11% (bkg) © 5% (lum) , 
Final state 3 (tjj): 1.5% (stat) © 11% (exp) © 9.1% (bkg) © 5% (lum) , (31) 

where in each final state the first term corresponds to the statistical error (including 
the background); the second is the experimental uncertainty from jet energy scaling, 
b tagging, etc.; the third one arises from background normalisation, and the last one 
from the luminosity determination. Errors are summed in quadrature, and amount 
to a 13.5%, 20.8% and 15.2%, respectively. If, in the real experiment, these errors 
can be decreased, e.g. with a jet energy calibration better than expected or with a 
different reconstruction method, then the results will be correspondingly better. The 
experimental uncertainty for the ratio R{t/t) in the final state 1 (tj) has been not 
estimated yet. We assume a 2% statistical error, consistent with the one in the first of 
Eqs. (l3Tl) . The luminosity uncertainty cancels in the ratio, which is also less sensitive to 
the background normalisation because the largest background is ti production, which 
contributes equally to t and i final states. Experimental errors from jet energy scaling, 
b tagging, etc. will likely affect the ratio to a lesser extent. For these uncertainties we 
tentatively assume a global value of 3%, to be summed in quadrature to the statistical 
error. Correlations between experimental systematic uncertainties in the different final 
states must be evaluated with a dedicated analysis, and have not been included in 
the limits. The correlation between a(t + i) and R(t/t) vanishes if the numbers of 
events observed N t , Nj follow Gaussian statistics with standard deviations y/N t , y/Nt- 
The statistical derivation of limits on anomalous couplings from observables is done 
in TopFit with the acceptance-rejection method, as explained in Refs. [5,6], and the 
limits presented here correspond to a confidence level of 68.3%. 

Before discussing the combined limits on the four Wtb couplings it is illustrative 
to consider examples in which some parameters are restricted to their SM values. 
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Figure 9: Estimated two-dimensional limits (with 68.3% CL) on (Vl, Vr) and (<7l, #.??), 
obtained from measurement of single top cross sections, with and without the ratio 
R{t/t) for the tj final state. 

In Fig. [9] we display two-dimensional limits in the case that all couplings except the 
pairs (Vl,Vr), (gL, go), respectively, take their SM values. We distinguish the limits 
obtained only from cross sections (in the three channels) and when the ratio R{t/t) is 
also included. The most interesting features of these results are: 

(i) Limits on {Vl, Vr) reveal that, as anticipated in section El the measurement of 
single top cross sections alone cannot discriminate Vl against Vr, and the ratio 
R(t/t) must be included as well. This ratio is equivalent to the asymmetry 
proposed in Ref. [2]. 

(ii) The inclusion of R(i/t) gives a moderate improvement in the limits on (gL,gn)- 

The combined limits obtained leaving the four parameters arbitrary are presented 
in Fig. [lOl Bounds on Vl and Vr are rather loose even including R(i/t), because of 
cancellations with terms involving §l and especially gR. Limits on gL and gR are also 
weaker than the corresponding ones in Fig. [9l up to a factor of two. 

Some comments regarding these limits are in order. As it can be clearly seen from 
the left plot in Fig. [TOl single top cross section measurements will not provide a mea- 
surement of Vl by themselves, unless additional assumptions on the rest of couplings 
are made. Notice that even setting Vr = 0, as it is often done in the literature, gives 
a large interval [0.5, 1.2] for Vl at 68.3% CL. The large allowed range for Vl is caused 
by the (partial) cancellations among the contributions to the cross sections involving 
Vl-> Vr, gL arid gR, and the limited experimental sensitivity. In this way, sets of these 
parameters very different from the SM values Vl — 1, Vr — gL — gR — give approxi- 
mate agreement with the SM predictions for cross sections in the three channels. The 
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Figure 10: Combined limits on Wtb couplings from single top cross section measure- 
ments, including R(t/t). The two graphs correspond to different projections of the 
4-dimensional allowed region (with 68.3% CL). 

weakness of the limits obtained emphasises the importance of a combination with top 
decay observables, which will be carried out in the next subsection. 

6.2 Combination with top decay observables 

Angular distributions and asymmetries in top decays are in many ways complementary 
to single top cross sections, and the combination of both provides much stronger limits 
on Wtb couplings. While the former cannot fix the value of Vl, they are much more 
sensitive to Vr, gi, and qr than the latter. We will not address here the interplay 
among single top production and top decay distributions in detail. This discussion 
will be presented elsewhere, when all the details regarding the experimental sensitivity 
for the observables involved will be studied. Still, it is very interesting to know the 
possible result of a combination, not only because of the better limits obtained, but 
also to test whether top decay asymmetries help improve the limits on Wtb couplings 
so that the efficiency variations in single top production due to anomalous couplings 
are small. 

In this combination we assume for the observables A± and p^i the experimental 
sensitivities of top pair production. These observables are redundant but they can be 
combined as long as their correlation matrix is not singular [6]. The combination of 
A±, Pr,l measurements in top pair and single top production will certainly improve 
the sensitivities with respect to these values and give tighter constraints on anomalous 
couplings. For the asymmetry ratio ru we tentatively assume an experimental error 
of 2%, uncorrelated with the former. Results are not strongly dependent on the ru 
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precision, as it will be demonstrated at the end of this subsection. We do not include 
r v \ in the fits, because it gives no improvement over alone unless their experimental 
uncertainties are similar. This is unlikely, due to the difficulty in reconstructing the 
momentum of the missing neutrino. For the combination we do not include R(t/t) 
either. This observable is superseeded by A± and Pr,l, and it would add some extra 
information to the fits only if its experimental precision was much better than the one 
estimated in the previous subsection. 

The combined limits obtained from the expected measurements of single top cross 
sections, A±, p R<L and are presented in Fig. HH These limits are far better than the 




a, A p,r 



Figure 11: Combined limits on Wtb couplings from single top cross section measure- 
ments (excluding R(t/t)) and top decay observables A±, pR t L, ru- The two graphs 
correspond to different projections of the 4-dimensional allowed region (with 68.3% 
CL). 

ones obtained either with single top cross sections or top decay observables alone, and 
clearly show the benefit of the combination among them. In particular, 

(i) Vl is bounded with a relatively good precision, 0.85 <Vl< 1.11 at one sigma, 
only a factor of 1.5 worse than the limit 0.92 < Vl < 1.10 that can be ob- 
tained from single top cross section measurements under the assumption that all 
anomalous couplings vanish. 

(ii) The constraints on Vr and gi are moderately strong, due to a fine-tuned can- 
cellation between them in A± and Pr ; l- This cancellation is decreased by the 
measurement of r w with a precision of ~ 8% or better, as argued below. 

(iii) Limits on gn are very stringent, —0.012 < g R < 0.024, even as good as the ones 
which have been previously obtained from top decay observables for Vl = 1 and 
assuming no cancellation between Vr and gi [6]. 



28 



We finally point out the important role of the asymmetry ratio r^, introduced in 
this paper, in order to improve limits on Wtb couplings. In Fig. [12] we display the 
projection on the (Vr, gi) plane of the 4-dimensional combined limits, distinguishing 
the cases when ru is not measured or it is measured with precisions of 8% and 2%. 
In the first case we can clearly observe the fine-tuned cancellation between Vr and gi 
present in the observables A± and Pr,l- This cancellation was already pointed out 
in Ref. [6]. Single top cross section measurements do not significantly reduce it, but 
when the ru measurement is added to the fits the cancellation is greatly decreased, 
even if the precision in the measurement is not very high. Thus, the limits presented 
in Fig. [TT] will likely be maintained to a large extent even if the 2% goal for r bt cannot 
be achieved. 




Figure 12: Projection on the (Vr^l) plane of the combined limits on Wtb couplings 
from single top cross section measurements and top decay observables A±, Pr,l, without 
r b i and with precisions of 8%, 2%. 



6.3 Efficiency variation with anomalous couplings 

With the limits in Figs. CCD and [[2] the efficiency variation in single top production can 
be evaluated, using some benchmark points. We select the following parameter sets, 
with anomalous couplings in the boundary of the 68.3% CL regions: 

Set A V L = 1, V R = 0.3, g L = 0.15, 

Set B V L = 1 , g R = 0.024 , (32) 

with the rest of couplings not explicitly written taken to zero. Events are generated 
for the tj, tbj, tb, tW~ and tW~b processes (and their charge conjugate) with each set 
of couplings, and simulated as before. The selection criteria (1-3) defined to isolate 
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t-channel, tb and tW~ production, respectively, are applied. The resulting efficiencies 
are gathered in Table QUI for the SM and the two sets of anomalous couplings. The 
quoted errors correspond to the Monte Carlo statistical error. For sets A and B effi- 

Final state 1 (tj) 

SM Set A Set B 

tj 2.02 ±0.02 2.08 ±0.02 2.07 ±0.02 

tbj 1.44 ±0.02 1.40 ±0.02 1.46 ± 0.02 

tb 0.76 ±0.01 0.73 ±0.01 0.74 ± 0.01 

tW 0.173 ±0.008 0.161 ±0.007 0.173 ±0.008 

tW'b 0.146 ±0.004 0.150 ±0.004 0.151 ±0.004 



Final state 2 (tb) 

SM Set A Set B 

tj 0.102 ±0.005 0.103 ±0.005 0.104 ±0.005 

tbj 0.328 ±0.009 0.315 ± 0.009 0.328 ± 0.009 

tb 2.70 ±0.03 2.75 ±0.03 2.70 ± 0.03 

tW~ 0.073 ±0.005 0.082 ±0.005 0.065 ± 0.005 

tW~b 0.071 ±0.003 0.071 ±0.003 0.071 ± 0.003 







Final state 3 (tjj) 




SM 


Set A 


Set B 


tj 


0.146 ±0.006 0.158 ±0.006 


0.132 ±0.006 


tbj 


0.294 ±0.00£ 


\ 0.313 ±0.009 


0.292 ±0.008 


tb 


0.33 ±0.10 


0.38 ±0.10 


0.32 ±0.10 


tw~ 


2.38 ±0.03 


2.47 ±0.03 


2.39 ±0.03 


tW~b 


2.37 ±0.02 


2.38 ±0.02 


2.36 ±0.02 



Table 10: Efficiencies (xlOO) for the different single top production processes in the 
three final states studied (after selection criteria), within the SM and for two sets of 
anomalous couplings. 

ciency variations are rather mild for the relevant contributions. (Processes contributing 
marginally have larger variations partly due to statistics, but these variations are ir- 
relevant for the limits finally obtained.) In the tj final state, t-channel processes have 
efficiency variations smaller than 3%, well below the 13.5% experimental uncertainty. 
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In the tb final state the tb and tbj processes have a variation up to 4.1%, which is also 
smaller than the experimental error of 20.8%, and in the tjj final state the tW~ and 
tW~b have a maximum variation of 3.8%, to be compared with the 15.2% experimental 
error. Hence, considering the efficiency as the SM one is a good first approximation, 
although the efficiency variation can be taken into account in a more detailed analysis. 

7 Conclusions 

Models beyond the SM allow for new physics effects in the Wtb vertex, either devia- 
tions from unity in the CKM matrix element V t b [50,51] or Wtb anomalous couplings, 
generated radiatively [52, 53] or from new physics at a higher scale. The main pur- 
pose of this paper has been to obtain expressions for single top cross sections at LHC 
involving arbitrary Wtb couplings. For definiteness we have worked within the frame- 
work of gauge invariant effective operators, which allows to relate the Wtb and gWtb 
couplings from gauge invariance. SM extensions at the electroweak scale can also give 
radiative corrections to both triple and quartic vertices, and in principle the relation 
between Wtb and gWtb may not be exactly the one predicted by the SU(2)l x U(l)y 
symmetry, which is broken at low energies. However, in this case the analysis becomes 
model-dependent and so the framework of gauge invariant effective operators remains 
simpler. 

We have explicitly shown that the only relevant Wtb couplings for single top produc- 
tion are the usual 7^ and a^ u q v terms, despite the fact that in some of the production 
processes the top and/or bottom quarks involved in the Wtb interaction are far from 
their mass shell. For this, we have introduced a general parameterisation of the Wtb 
vertex for off-shell top and bottom quarks, in terms of the on-shell Lagrangian (with 
7 M and (T^ u q u terms) plus two "off-shell" operators Oi 2 which vanish when the top and 
bottom quarks are on their mass shell. Assuming that these couplings arise from gauge 
invariant operators, have quartic gWtb vertices O^fl associated by gauge symme- 
try, which also contribute to single top production. We have seen that the combined 
contributions involving and identically cancel, making the single top cross 
sections independent of the off-shell couplings. This cancellation reduces the number 
of relevant parameters in the Wtb vertex from 8 to 4 and thus simplifies setting lim- 
its from single top cross sections. Furthermore, it also implies that Gordon identities 
can actually be used on gauge invariant new physics contributions to the Wtb vertex, 
rewritting and a^ u k^ terms even for off-shell t and b. 

We have provided expressions for single top cross sections in terms of Wtb anoma- 
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lous couplings for (i) the t-channel process, including tj and tbj, for which a matching 
has been performed to avoid double counting; (ii) s-channel tb production; (iii) tW~ 
associated production including the tW~b correction, from which the resonant ti con- 
tribution is subtracted at the cross section level. These cross sections have been written 
as a sum of products of Wtb couplings times numerical coefficients, 

o = a S M(V? + K V «V* + K V L v xV L V R + K^gl + K 9 *g 2 R + . ..) , (33) 

with the couplings defined in Eqs. (pQ) and the dots standing for further interference 
terms. These expressions can be used to fit the LHC measurement of single top cross 
sections and simultaneously obtain a measurement of V t b and bounds on anomalous 
couplings. We have shown how this could be done, taking into account the theoretical 
uncertainty in single top cross sections and k factors, as well as the expected experi- 
mental uncertainty in cross sections [47]. A difference with respect to usual single top 
analyses is that here all single top processes are considered as signals, with different 
efficiencies and also with a different dependence on anomalous couplings. Indeed, this 
seems the most reasonable approach in an analysis aiming to measure Vtb and anoma- 
lous couplings, since the cross section of all single top processes depend on them. 
Backgrounds are then constituted by processes with cross sections independent of the 
Wtb interaction, as for example ti and W/Z production plus jets. 

The limits obtained using only single top cross sections are not very strong because 
of the experimental uncertainties assumed, of 13%, 21% and 15% for the t-channel, tb 
and tW~ processes, respectively (which may be a little conservative). Moreover, can- 
cellations among the different couplings, which have not been considered in previous 
literature, prevent us from obtaining precise bounds using only cross section measure- 
ments. Nevertheless, limits can be greatly improved with the combination with top 
decay observables, such as the angular asymmetries A±, helicity ratios pr^l [5] and a 
spin asymmetry ratio ru introduced here. We have performed the combination using 
the expected precision for A± and pr^ in top pair production [6] and an estimate of 
the sensitivity in the ru measurement. It must be remarked that the limits obtained 
make no assumption on the Wtb couplings (except that they are CP-conserving). Due 
to the non-linearity of the equations used and the limited precision of the experimen- 
tal observables, determining the Wtb couplings without ambiguities is more involved 
than merely counting parameters and finding an equal number of independent observ- 
ables. Actually, it is non-trivial to find stringent limits on Wtb couplings avoiding the 
fine-tuned cancellations that can occur, mainly between Vr and Ql- The asymmetry 
ratio ru (or any equivalent observable) plays a key role in improving the limits, as 
demonstrated in Fig. [T2l 
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The results obtained from the global fit are very promising. Precise bounds on 
several couplings can be achieved, as well as a measurement of V tb which is only a 
factor of 1.5 less precise than the (model-dependent) one obtained setting all anomalous 
couplings to zero. These results have been shown to depend weakly on the sensitivity 
to ru if it is better than ~ 8%, which is likely to happen given the sensitivities of other 
top decay observables. This combination deserves a detailed investigation when the 
expected experimental precision of all the observables involved is known in one or more 
top production and decay processes, and it will be presented elsewhere. 
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A Cancellation of off-shell operator contributions 

In this appendix we explicitly prove the cancellation among off-shell triple and quartic 
contributions to the gb — > tW~ and gg — > tW~b amplitudes. We first rewrite in 
Eqs. (fTTj) expanding the a^ u terms in both operators, obtaining 

Oi = -- / J^ r b[-(f 1L P L + f 1R P R )Y^t-m t ) 

+ (^-m b )Y(fiLP L + fmP R )} t W~ + H.c. , 

+ {^-m b )Y{f 2L P L + f m P R )} t W~ + H.c. , (34) 

where p t and p b are the momenta of the quarks involved in the vertex, following the 
fermion flow. Note that in deriving Eqs. ( 1341 ) from the definition in Eqs. ( 1X71 ) we have 
not used Gordon identities, so that the two sets of equations are completely equivalent 
for t, b off-shell. These alternative expressions for 0± i2 are extremely useful to prove 
the cancellation among diagrams. Incidentally, Eqs. fl34l make apparent the fact that 
01,2 cancel when both the top and bottom quarks are on their mass shell, and they 
also show that these two operators are not independent if either the top or the bottom 
quark is on-shell. For top on-shell, 

fin) - 2 (f 2L = f 1L , f 2R = f 1R ) = , (35) 



33 



whereas for the bottom quark on-shell 

fa) + 2 (f 2L = f 1L , f 2R = f 1R ) = , (36) 

in obvious notation. Here we will restrict ourselves to the cancellation among 2 and 
02 for brevity, because the proof for 0\ and is almost identical, as it can be 
observed from their expressions in Eqs. ( 1341) . 



A.l Cancellation in the gb — > tW amplitude 

There are three Feynman diagrams contributing to this process, the two ones in Fig. [5], 
involving a Wtb vertex, and a third one with a quartic gWtb coupling. We denote by 
Pi? P2, P3 and p4 the momenta of the external gluon, b quark, top quark and W boson, 
respectively. The matrix element for the t-channel diagram reads 

Mx = -^-^^u(P3)l U (^ + rn t )m-m t )r(f2RPL + f2LP R )}u(p 2 ) 

V 2Mw Pb ~ m t 

xe*(p 4 )e„(Pi), (37) 

using the expression for 2 in Eqs. (i34l ) and the fact that the b quark is on-shell, 
p/ 2 u(p 2 ) = rrii } u(p 2 ). The momentum of the internal top quark is p§ = p 2 — p^. The 
factor between brackets corresponds to the 2 vertex, up to numerical constants. We 
have suppressed colour indices and a A a /2 factor, which are common to the three 
diagrams, and after colour averaging and summing amount to a global factor of 1/6. 
The internal top quark propagator cancels with the (pi —mt) vertex factor and the 
matrix element simplifies to 

Mx = -^-u(p 3 )YY(f2 R PL + f2LP R )u(p 2 )e;(p 4 )e l/ (p 1 ). (38) 
The matrix element for the s-channel diagram reads 

M 2 = ^^^^«(p 3 )[(/ 2fl P i + / 2L P fi ) 7 M(y 6 -m 6 )](y 6 +m fe )7^(p 2 ) 

V 2Mw Pe ~ m b 

x<(P4)e„(Pi), (39) 

where p& = pi + p 2 is the momentum of the internal b quark, and we have used that 
the external top quark is on-shell, uip^p 1 ?, = u{pz)m t . This simplifies to 

M 2 = -^-u(p 3 )Yl u (f2RPL + f2LP R )u(p 2 )e;(p 4 )e l/ (p 1 ). (40) 
Finally, the diagram involving the quartic interaction gives 

M( 4) = _^^^ 3) ^ (/2R p L + /2L p fl)M(p2)e * (p4) ^ (pi)) (41) 
IVlw 

and the cancellation follows from {7 M ,7 I/ } = 2g fiu . 
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A. 2 Cancellation in the gg — > tW b amplitude 



The diagrams with Wtb vertices contributing to this process are the same as in the SM, 
shown in Fig. [6l and the additional diagrams with quartic gWtb vertices are displayed 
in Fig. HH We number the diagrams from left to right and top to bottom. There are 
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Figure 13: Feynman diagrams for gg — > tW b involving anomalous quartic vertices. 
The two initial state gluons are labelled as g± and gi for clarity. 



two different colour flows contributing to the amplitude. The two diagrams involving 
a triple gluon vertex contribute to both colour flows but cancel with the last diagram 
in Fig. HH 

M 5 + Ms + Mf = , (42) 



in the same way as for gb — > tW~ , explained in the previous subsection. Then, we 
concentrate ourselves on the remaining diagrams. The momenta of the two initial state 
gluons are p\ and p%, and the momenta of the top quark, W boson and b quark are p 3 , p 4 
and ps, respectively. Inserting the expression of O2 from Eqs. (I34L using the equations 
of motion for the external top and bottom quarks and simplifying the propagators, the 
matrix elements for the remaining diagrams are 
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+ 2 1 2 U( P3 ) (f 2R P L + f2LP R h' 7 l' 1 (]/lO+m b )Y V(p 5 ) 
PlQ ~ m b 

M 3 = -J^-^^u(p 3 ) (f 2R P L + h L P R )n v (l/7+m b h a v(p 5 ) x e, 

V 2Mw Pi — rn b 

M 4 = 2 1 2 u( P3 ) (f 2R P L + f 2L P R )YY(l/xo+m b )Yv(p 5 ) x e , 

V 2M W Pio - m b 

y ZM W p 9 — m t 

M 7 = -Utt- «(P3) 7 i/ (y6+m t )7V(/2i ? P i + «(p 6 ) x e , (43) 



^ pi - 



where e stands for e* (p4)e J/ (p 2 )ecr(Pi)- The additional momenta introduced are p 6 
Pa—P2,Pr = Pi—P5,P9 = Pa-Pi, Pio = P2 ~ Pb- The first four diagrams in Fig. 
give 

M? = -^p -^-^u(p 3 ) Y(^+m t )g^(f 2R P L + f 2L P R ) v(p 5 ) x e , 

% £>6 — m t 

M W = _V2^gl r{]/9+mt)g ^ (fm p L + /2L p fl) u(p5 ) x e ; 

iW w P9 — m t 

Mf = -j-L^fifo) fePL + / 2 LPi ? )^(y 7 +m 6 ) 7 <T «(p B ) x 6 , 

iw w p 7 — m b 

M (A) = _V2ggl 1 (/ ^ Pl + / 2L p i? )^( 2 / 10+m6 ) T ^( P5 ) x e . (44) 

Using {7^, 7^} = 2g Ml/ , for one colour flow we have 

M! + M 3 + M 7 + m¥ ] + M { 3 4) = 0, (45) 

and for the other one 

M 2 + M 4 + M 6 + Mf + A4f = . (46) 

Finally, we point out that introducing a gauge non-invariant regulator for the top 
quark pole singularities, like the usual prescription for the propagators p 2 — m 2 — ► 
p 2 — m 2 + im t T t , spoils the cancellation. Alternative gauge-invariant prescriptions can 
be found for example in Ref. [45]. 



B Top decay observables 

For completeness, we briefly introduce here the definitions of top decay observables used 
in the combination with single top cross sections. An extended discussion regarding 
these observables can be found in Ref. [5]. 
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The charged lepton angular distribution in the W rest frame is very sensitive to Wtb 
anomalous couplings. It is defined in terms of the angle 6} between the charged lepton 
in the W rest frame and the W momentum in the top quark rest frame. An important 
point is that this distribution is independent of the top quark production mechanism, 
in particular independent of the top polarisation, and thus the related observables can 
be defined and measured either in top pair or single top production. The normalised 
angular distribution of the charged lepton can be written as 

l Yi^¥rl [l+me ' fFa+ l (l - m e " ,? Fl + 1 sm2 e "' Fa ■ m 

The helicity fractions Fj = Ti/Y are the normalised partial widths for the top decay to 
the three W helicity states (we drop here the t subindex in the top width for simplicity). 
A fit to the cos 9} distribution allows to extract from experiment the values of Fj, which 
are not independent but satisfy Fr + Fl + Fq — 1 by definition. From the measurement 
of helicity fractions one can constrain Wtb anomalous couplings [29]. Alternatively, 
from this distribution one can measure the helicity ratios 

_ F RL F RL 

PR,L = -p- = -£T- , (48) 
J- -TO 

which are independent quantities. For any fixed z in the interval [—1, 1], one can also 
define an asymmetry 

= N(cos6*>z)-N(cos6*<z) 
z N(cos6* > z) + iV(cos0; < z) ' { ' 

The most obvious choice is z = 0, giving the forward-backward asymmetry Afb [32,37]. 
But more convenient choices are z = =f(2 2 / 3 — 1) [5]. Defining (3 = 2 1 / 3 — 1, we have 

2 = _(2 2 /3_ 1) ^ A g = A + = 3/3[F + (l + /3)F R ], 
z = (fl^-\) _ A Z = A^ = -3P[F + (1 + P)F L ]. (50) 

These asymmetries are obviously determined by the W helicity fractions (or ratios) 
and conversely, from their measurement the helicity fractions and ratios can be recon- 
structed. Moreover, the charged lepton energy distribution in the top quark rest frame 
is uniquely determined by the helicity fractions and the top, bottom and W boson 
masses. Among these observables we select A± and Pr^l, which are the most sensi- 
tive to Wtb anomalous couplings, for the combination with single top cross sections. 
They take the tree-level SM values A + = 0.5482, A_ = -0.8397, p R = 5.1 x 10" 4 and 
Pl = 0.423, for a top quark on its mass shell. 

Further observables can be built involving the top spin. For the decay t — > W + b — > 
£ + ub,qq'b, the angular distributions of X = £ + ,i / ,q,q',W + ,b (which are called "spin 
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analysers") in the top quark rest frame are given by [30,31] 

1 dT 1 



(1 + ax cos 6x) (51) 



T d cos 9x 2 

with 9x the angle between the three-momentum of X in the t rest frame and the top 
spin direction. The constants ax are called "spin analysing power" of X and can range 
between —1 and 1. For the decay of a top antiquark the distributions are the same, 
with a x = —oix as long as CP is conserved in the decay. In the SM, a^+ — aq' = 1, 
a v = a q = —0.319, otw+ = — Oi b = 0.406 at the tree level (q and q' are the up- and 
down-type quarks, respectively, resulting from the W decay). Tree-level expressions of 
the spin analysing power constants for a CP-conserving Wtb vertex with t, b on-shell 
as in Eq. (jTJ) have been obtained in Ref. [5] within the narrow width approximation. 

In the t-channel single top process the top quarks are produced with a high de- 
gree of polarisation along the direction of the final state jet [4]. The corresponding 
distributions are 

-{1 + Pa x cos6 x ), (52) 



TdcosO x 2 

where the angles 9x are now measured using as spin direction the jet three-momentum 
in the top quark rest frame, and P ~ 0.89 is the top polarisation along this axis. 
(This value, calculated with the matching of tj/tbj at 10 GeV, is in good agreement 
with the polarisation calculated with the subtraction method [4].) Forward-backward 
asymmetries can be built using these distributions, 

= N (cos Ox > 0) - N (cos Ox >0) = 1 , . 

X N(cos0 x > 0) + N(cos0 x > 0) 2 &X ' { ' 

These spin asymmetries depend on the top polarisation P but their ratios do not, for 
example 

A b a b 

-r = — = r bi , 
A e a e 

Ai a e 

only depend on anomalous couplings and the masses involved. Additionally, their ex- 
perimental systematic errors are expected to be smaller than for the spin asymmetries 
themselves. As it happens with A± and Pr,l, these asymmetry ratios can be mea- 
sured in top pair production as well. In tt production the top quarks are unpolarised 
at the tree level, but their spins are correlated [54,55]. Suitable ratios of spin cor- 
relation asymmetries can provide a measurement of r b i and r v \ but the precision will 
likely be worse than in single top production because spin correlation asymmetries are 
numerically smaller, e.g. by a factor of 8 in the tt semileptonic decay channel. 
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C Effect of top quark decay 



Here we examine the effect of top decay in the k factors, concentrating on the t- 
channel process for brevity. We calculate the k factors for the processes with top decay 
numerically, following the procedure explained in section [5] for tW~b production. For 
the processes with top on-shell the calculation is carried out as in section [3l with the 
difference that here for the sake of computational speed we perform the calculations 
for tj without the veto on p^ SR , which is sufficient to illustrate the differences. Results 
for the quadratic terms are shown in Table HU As expected, the effect of the top 
width is small, at the 1 — 2% level, and in most cases it is smaller than the theoretical 
uncertainty. The differences have little numerical relevance in the limits obtained for 
anomalous couplings. 

tj tbj tj ibj 

on off on off on off on off 

V% 0.895 0.888 0.927 0.917 1.101 1.093 1.068 1.064 

g\ 1.47 1.47 1.96 1.96 2.12 2.13 2.98 3.01 

g% 2.00 1.98 2.97 2.95 1.63 1.62 2.08 2.07 

Table 11: k factors for the quadratic terms in the tj, tbj processes and their charge 
conjugate, calculated for top on-shell or decaying t — > W + b — ► tvb (labelled as "on" 
and "off", respectively). 
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